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ABSTRACT 

Context. It has been shown that the luminosity of active galactic nuclei and the size of their broad line region obey a simple relation of 
the type Rblr = aL''', from faint Seyfert nuclei to bright quasars, allowing single-epoch determination of the central black hole mass 
Mbh = bUAjj^ from their luminosity L and width of H/j emission line. Adopting this mass determination for cosmological studies 
requires the extrapolation to high redshift and luminosity of a relation whose calibration, relies so far on reverberation mapping mea- 
surements performed for L S10''*erg s"' and redshift z ~ 0.4. 

Aims. We initiated a campaign for the spectrophotometric monitoring of a few luminous, intermediate redshift quasars whose apparent 
magnitude, V < 15.7, allows observations with a 1.8 m telescope, aimed at proving that emission lines vary and respond to continuum 
variations even for luminosities ^10**^ erg s"', and determining eventually their M^h from reverberation mapping. 
Methods. We have repeatedly performed simultaneous spectrophotometric observations of quasars and reference stars to determine 
relative variability of continuum and emission lines. We describe the observations and methods of analysis. 

Results. For the quasars PG 1634+706 and PG 1247+268 we obtain hght-curves respectively for CI1I](^^1909A), MglI(i^2798A) 
and for ClV(/i/il549A), C1I1](/1/11909A) emission lines with the relevant continua. During 3.2 years of observation, in the former case 
no continuum variability has been detected and the evidence for line variability is marginal, while in the latter case both continuum 
and line variability are detected with high significance and the line variations appear correlated with continuum variations. 
Conclusions. The detection of the emission line variability in a quasar with L ~ 10*^ erg s"' encourages the prosecution of the moni- 
toring campaign which should provide a black hole mass estimate in other 5-6 years, constraining the mass-luminosity relation in a 
poorly explored range of luminosity. 

Key words. galaxies:active - quasars: emission lines - quasars: general - quasar:individual:PG 1634+706, PG 1247+268 



1. Introduction 

Supermassive black holes (SMBHs) are believed to inhabit most, 
if not all, the bulges of present-epoch galaxies (Kormendy & 
|Richstone] 1995| l, and strong evidences exist of a correlation be- 



tween the 



an accretion disk surrounding the black hole. Emission lines re- 
spond to variation of the ionizing continuum. Although the phys- 
ical origin of these variations is poorly known (e.g. Trevese & 
Vagnetti||2002||Vanden Berk et al.||2004||De Vries et al.||2005| l 



5lack hole mass and either the mass Mtmige and lu 



minosity (Marconi et al. 2003 and refs. therein) or the veloc- 



ity dispersion cr., of t he host bulge ( |Ferrarese & Merrit] |2000[ 



[Tremaine et al. 2002[ l. This strongly suggests that the formation 
and growth of SMBHs and galaxies are physically related pro- 
cesses and provides a basis for a theory of cosmic structure for- 
mation, incorporating the feedback from Active Galactic Nuclei 
(AGNs) (Sil k & Rees] [19981 |Vittorini, Shankar, & CavaEerel 
[2005 and refs. therein). Black hole masses determinations based 
on stellar or gas kinematics are intrinsically limited, by angu- 
lar resolution, to relatively nearby objects and cannot be applied 
to bright AGNs where the nuclear light prevails over the galac- 
tic component, just in the central region where the galactic gas 
or star motion is dominated by the black hole gravitational field. 
The reverberation-mapping technique does not suffer of this lim- 
itation and represents the only mean to measure the mass of 
SMBH in bright AGNs. 

Emission lines, in the optical-UV region, are interpreted as 
recombination of a gas which is photoionized by the continuum 
radiation emitted by the inner region of the nucleus, presumably 



it is possible to use the response of lines to continuum varia- 
tions to investigate the structure of the line emitting region. This 
requires long campaigns of accurate spectrophotomeiic monitor- 
ing of AGNs, which have led in the past to major progresses to- 
wards understanding the physics of the "atmosphere" of Seyfert 
1 galaxies. A summary of these results is given in |Peterson| 
( |TW3] l. 

Line widths, e.g. the FWHM A//^ of the emission line, 
coiTespond to r.m.s. velocities of the emitting gas clouds. A 
cross-correlation analysis of continuum and emission-line light- 
curves, evidencing a time delay t of line respect to continuum 
variations, allows to estimate the size R - t/c of the region 
where the line photons are generated. If the gas motion in the 
Broad Line Region (BLR) is dominated by gravitation (Peterson] 
& Wandel, ,2000) , the size estimate Rblr can be combined with 
the line width to yield a primary estimate of the virial mass of 
the black hole Mgn ^Ii^/GRblr and the relevant Eddington 
ratio. For Seyfert 1 galaxies typical BLR sizes are of the order of 
light-days to light- weeks. Similar studies are more difficult for 
quasars (QSOs) which require a longer monitoring. 
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A long term campaign for a subsample of 28 QSOs was 
started in 1991 with the Wise 1.0m and the Steward 2.3m tele- 
scope s ()Maoz et al. [1994 1. As a result, nine years later Kaspi 
|et al.| ( |2000| l provided mass estimates for the entire sample. The 
new data, combined with previous results on Seyfert 1 galax- 
ies, thus spanning a much wider range of intrinsic luminosity, 
allowed to establish an average relation between the intrinsic lu- 
minosity and the size, Rblr - ciU, with y ^ 0.7, which allows a 
secondary estimate of the black hole mass based on single-epoch 
observations of luminosity and line width: Mbh - bUA'^^^, 
where both the constant b and y are determined statistically from 
the available echo-mapping data. Recent studies show that y is 
in the range 0.5-0.7, depending on how luminosity is defined, 
which lines are selected for the echo-mapping and the fitting pro- 
cedure adopted ( Kaspi et al. 2005 Bentzetal. 2006jl. 

The extreme importance of secondary mass estimates relies 
on the fact that on the basis of single epoch observations it is 
possible to study the evolution in cosmic time of the mass distri- 
bution of QSOs/AGNs, and to extend the studies of the relation 
existing between QSOs and host bulges properties. However, 
the above correlations with primary masses, based on echo- 
mapping, were established for relatively close and faint AGNs 
with z < 0.4 and [^L^(5100A) ^ 10*^ erg s"'], thus it is 
presently unknown whether they can be extrapolated to higher 
luminosities and/or redshifts. For instance, the extrapolation of 
the Mbh - L relation (Kaspi et al. 2000| l, together with the 
assumption that the known Mbh - Mi,„ige - cr» relations holds 
( Tremaine et al.| [2002) , leads to predict the existence of galax- 
ies with Mbuige ~ 10'^ ' - 10''' '*Mq and cr, exceeding 800 km 
s Such galaxies have never been observed, and their exis- 
tence would put important constraints on galaxy formation mod- 



Table 1. The quasars monitored 



els (Netzer 2003 1. Therefore it is essential to extend the primary 
mass measures to higher redshifts and luminosities. On the other 
hand, for high QSO luminosities a large size of the broad line 
region is expected. This would cause both a smoothing of the 
line light-curve and larger time delays with respect to contin- 



uum variations ( Wilhite et al. 2005 1, thus the very detectability 



and the amplitude of line variations are open questions. 

A sample of objects with redshifts in the range 2 < z < 3.4 
and apparent magnitude as faint as niy ~18 is being monitored 
by Kaspi et al. ([2004 ) with the 9m Hobby-Eberly Telescope 
(HET; Ramsey et al. T1998 ) and new results have been presented 
recently ( ,Kaspi et al. ]2006^ . During their 6-year monitoring of 6 
QSOs, significant continuum and emission-line variations were 
detected in all targets and a preliminary black hole mass estima- 
tion is given for one of them. 

In the present paper we describe a new monitoring campaign 
limited to objects with V < 15.7 and 1 < z < 4 which, thanks 
to their apparent brightness, can be observed with the medium- 
small 1 .82 m Copernicus Telescope at Cima Ekar (Asiago, Italy), 
through a service mode scheduling of a long term monitoring, 
and allow to investigate whether: i) echo-mapping is feasible 
for objects as bright as /IL^(5100A) ~ 10"*^ erg s"' and ii) the 
^BL«-luminosity correlation can be extrapolated to such bright- 
ness. The paper is organized as follows. Section 2 describes the 
sample, observations and the data reduction procedure. Section 
3 describes the results for two quasars of the sample. Section 4 
summarizes the results and discuss future prospects. In the fol- 
lowing we derive /IL^(5100A) from the flux in the Johnson V 
band, extrapolating the flux density to the rest-frame A = 5100A 



with a power law fy 



-0.5 



, and assuming a standard cosmology 
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Fig. 1. Average spectra of PG 1634+706 (upper panel) and PG 
1247+268 (lower panel). Spectral ranges for continuum deter- 
mination (short ticks) and line flux evaluation (long ticks), as 
reported in Table 2, are shown. 



2. Observations 

2.1. Object selection 



The sample has been extracted from the 'Veron-Cetty & Veron 
(2003) (11th ed.) catalog with the condition 6 > 0, V < 15.7 
mag and z > 1 in order to select objects of bright enough intrinsic 
luminosity to investi gate the bright end extens ion of the Rblr vs. 



/iLj(5100A) relation ( Kaspi et al. 



2000 



2005 1. These conditions 



Ho = 70 km S-' Mpc-', = 0.3, and Qa = 0.7. 



identify 12 objects, only four of which were monitored, owing 
to the limits on observing time. These four objects are listed in 
Table 1. 

Observations at intermediate redshift allow to sample the 
variabihty of Mgll ^2798, CIII] ^1909, CIV ^1559 lines, instead 
of Ha, Hp, Hy observed in the low redshift study of Kas pi et al.| 
(2000 1. This allows to study either BLR at smaller distance from 
the center, or regions of the same sizes of those producing the 
Balmer lines, but using lines which respond to different part of 
the continuum spectrum. The main emission lines falling within 
the observed wavelength interval are indicated for two objects 
in Table 2. Figure [T] shows the average spectra of PG 1634+706 
and PG 1247+268. 

With respect to other QSO monitoring programmes, ours 
is the first which includes the Mgll line in part of the ob- 
served sources. This line is particularly significant, because i) 
past Seyfert 1 monitoring campaigns conducted with lUE have 
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shown its lag to be similar to that of Hp (Clavel et al. 1991 
[Reichert et al.|[T994l l; ii) this line is most often used to derive es 
timates of black-hole mass from single-epoch spectra of high-z 
QSOs ( [McLure & Jarvisj [2002[ ), because its width is tightly cor- 
related to that of Hp. Deriving a lag for the Mgll line would 
therefore allow to estimate the black-hole mass most consis- 
tently with respect to the results of Seyferts, which are mostly 
based on the monitoring of Hp. 

APM 08279+5255 is one of the most luminous known QSOs 
if its emission is considered isotropic. However it has been 
shown to be lensed by a foreground galaxy (Irwin et al. 1998 1. 
Three components, separated by a few tenths of an arcsec have 
been detected in near-infrared images obtained with Keck tele- 
scope and different models of the lensing field predict a few 
days delay between photometric variations of the components 



(Ledoux et al. 1998 



Ibata et al. 1999 Egami et al. 2000 1: 



a short time compared with the expected time scale of intrin- 
sic variations. APM 08279-1-5255 is also a Broad Absorption 
Line QSO (I rwin et al.| |1998| l, which makes more difficult to 
define regions free from either emission or absorption features 
to measure continuum variations. The analysis of this object is 
deferred to a forthcoming paper HS 2154+2228 has been ob- 
served so far only 5 times and the analysis requires further mon- 
itoring. The other two objects PG 1247+268 and PG1634+706 
are analyzed in this paper with the aim of : i) verifying the ad- 
equacy of the observational data and reduction procedures, un- 
der the assumption that variability amplitudes and characteristic 
time scales can be extrapolated from the properties of fainter 
objects; ii) possibly detecting line variations in objects as lumi- 
nous as /lL^i(5100A) ~ lO**^ erg s"', and compare their amplitude 
with continuum variations. For both objects a star of comparable 
magnitude, as close as possible to the QSO, has been selected 
for the relative spectrophotometric calibration, described in the 
next section. 



2.2. Spectrophotometric Observations and Data Reduction 



Observations were carried out at Asiago 1.82 m telescope 
equipped with the Faint Object Spectrograph & Camera AFOSC 
which is a focal reducer with reduction factor of 0.58, designed 
to allow a quick switching between spectroscopic and imaging 
modes. The scale at the focal plane is 21.7 "/mm. The detector 
is a 1024x1024 thinned CCD array TEK1024 with 22x22 jum^ 
pixels corresponding to a scale of 0.473 arcsec pixel"' and a 
FOV of 8.14x8.14 arcmin^. We adopted a 8".44-wide sHt and 
a grism with a dispersion of 4.99 A pixel providing a typ- 
ical resolution of ~ 15 A in the spectral range 3500-8450 A. 
Spectrophotometric exposures are performed after orienting the 
slit to include both the QSO and the reference star of comparable 
magnitude, located at (12:50:11.5 +26:33:32) and (16:34:57.4 
+70:32:49) (J2000) forPG 1247+268 and PG 1647+706 respec- 
tively. The reference stars are included as internal calibrators for 
QSO spectra, as described by Maoz et aT] ( |1990 ) and Netzer et 
[HI] ( |1990) l. The wide slit is necessary to avoid different fractional 
losses of the QSO and star light due to possible non perfect slit 
alignment and differential refraction, which could cause spuri- 
ous variation of the flux ratios. Lamp flats and Hg-Cd ark spectra 
were also taken for wavelength calibration. At each epoch, typi- 
cal science observations consist of two consecutive exposures of 
1800 s. 



QSO and star spectra,2(/l) and S(A) were extracted with the 
standard IRAF procedures. The QSO/star ratio as a function of 
wavelength is computed for each exposure k - \,2 



(1) 



This quantity is independent of extinction changes during the 
night. This allows a check of consistency between the two expo- 
sures and the rejection of the data if inconsistencies occur (under 
the assumption that QSO variations are negligible on 1 h time 
scale). In fact, typical spectra of two consecutive exposures have 
a ratio fi^^\A)/f/^\A) of order unity, with deviations smaller than 
0.02 when averaged over 500 A, at least in flie 4000A - 7000A 
range. When discrepancies are larger than 0.04 both exposures 
are rejected. 

It is important to notice that, because of the large slit width, 
there are small changes in the A scale due to changes of the ob- 
ject position within the slit (which are in general negligible in 
the case of pairs of consecutive exposures). This suggested to 
proceed as follows. As a first step we register the zero point of 
the A scale. However, since the changes of position within the slit 
do not correspond exactly to rigid shifts of the A scale, we select 
portions of the spectra in wavelength intervals \AA\ :$ lOOOA 
around each of the QSO emission lines considered (see Table 2), 
and we determine the shifts for each portion. 

Once the wavelength scales are registered, the QSO and star 
spectra taken in the two exposures are co-added and the ratio 



|(2)^ 



(2) 



is computed, at each epoch /. The reference star is flux calibrated 
at a reference epoch. Since our aim is to compute relative flux 
variations, we reduce all quasar spectra to this reference epoch, 
by multiplying all the i-ii(A) by the flux calibrated star spectrum 
f^(A): f^(A) = jii(A)f^(A). We stress that we are not interested 
in the absolute flux calibration, whose accuracy is of the order 
of 20% and which is applied for the sole purpose reporting the 
QSO spectra in physical units. The star spectrum /^(A) adopted 
is the same for all epochs, thus it does not affect relative flux 
changes we want to measure. 

The QSO spectrum f^{A) - c{A) + 1{A) can be decom- 
posed in a line 1(A) and continuum c{A) spectra. Two values 
of the continuum, Csimn and cio„g at shorter and longer wave- 
length with respect to the most prominent QSO emission lines, 
are evaluated in regions as free as possible of other emission 
features, defined by the wavelengths ranges {Ashort,i,Ashori,2) and 
{Aiong,i,Aio„g;2) respectively. Table 2 reports the observer-frame 
wavelengths defining the continuum regions and the intervals for 
line integration. These were chosen on the basis of the analysis 
of [Clavel et aL] ( [T99T] l; |Reichert et aL] ( [T994) of UV spectra of 
low redshift AGNs observed with the International Ultraviolet 
Explorer (lUE), with slight modifications to maximize the S/N 
ratio {ciong of MgII(/l2798 A) falls outside the wavelength range 
covered by the lUE spectrograph). 

Line fluxes are computed as: 



//= { '[F''Q\A)-c'"'(A)]dA, 
Ja, 



(3) 



where c""(A) is the linear interpolation through c short and c;„„g. 
The extremes Ai and A2, also listed in Table 2, not necessarily 
coincide with Ashon = {Ashort,\ + A,,h„rt,2)/2 and Ai„„g = (Aio„g,i + 
Aio„g 2)l2, and are chosen to optimize the // signal to noise ratio. 

Direct B and R images in 15x15 arcmin fields centered on the 
QSOs were taken to check possible variability of the reference 
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Table 2. Wavelength intervals for lines and continua [A] 



Object 




line 


^ short, 1 


'^slwrl.2 


^1 


A2 






PG1634+706 


1.337 


C///] .11 909 


4252 


4298 


4382 


4590 


4630 


4674 


PG1634+706 


1.337 


Mg/W2798 


6215 


6285 


6251 


6777 


6755 


6815 


PG1247+268 


2.042 


C/V^1549 


4380 


4502 


4598 


4801 


5202 


5262 


PG1247+268 


2.042 


C///]^1909 


5535 


5595 


5770 


5900 


6025 


6085 



Clll] 



Mgll 



Clll] 



2800 3200 3600 
t (JD-2450000) 



2800 3200 3600 
t (JD-2450000) 




2800 3200 3600 
t (JD-2450000) 



2800 3200 3600 4000 
t (JD-2450000) 



Fig. 2. Light-curves of PG 1634+706 in the observer's frame. Fig. 3. Light-curves of PG 1247+268 in the observer's frame, 

as relative flux variations A///. Upper panels: emission lines, as relative flux variations A///. Upper panels: emission lines, 

middle panel: shorter wavelength continuum Cshon, lower panel: middle panel: shorter wavelength continuum Cshon, lower panel: 

longer wavelength continuum cio„g, left: CIII]//11909 A, right: longer wavelength continuum c/o„g, left: CIV//11549 A, right: 



MgH A2798 A. 



CIII] .11909 A. 



star. Comparison of the reference star with the brightest star in 
the field showed a r.m.s. fractional flux variation of 0.02 in both 
cases. Thus we can assume that the reference star is not variable, 
at this level of accuracy. 

Relative-flux calibration 1-cr errors, reported in figures|2]and 
[3] are estimated by adding in quadrature the r.m.s. uncertainty 
on the flux of the reference star to the r.m.s. fractional uncertain- 
ties of line or continuum fluxes, which, in turn, are computed by 
comparing pairs of exposures which are added to form quasar 
spectra of individual epochs. 



3. Results and discussion 

Figures|2]and|3]show the light-curves in the observer's frame, for 
emission lines and continua listed in Table 2, for PG 1634+706 
and PG1247+268 respectively. The light curves are expressed 
in terms of the relative flux variations A///, with respect to the 
flux / at the reference epoch where the absolute calibration was 
performed (see Section 2.2). The total time base is 3.3 years in 
the observer's frame, corresponding in the rest-frame to 1.4 yr 
for the former QSOs and 1.1 yr for the latter. In the case of 
PG 1634+706 no significant variations are detected: neither in 
CIII] line, which has the lowest S/N ratio due to its intrinsic 
faintness, nor in the Mgll line. The relevant continua are also 



consistent with no variability at a level of A/// ~ 0.02-0.03 
r.m.s.. 

The result is different in the case of PG 1247+268. The 
faintest line CIII] /11909 A shows a marginal evidence of vari- 
ability and the stronger line CIV /1 1549 A seems to decrease 
steadily during the observing period. 

An almost steady decrease does not allow to derive possi- 
ble line-continuum time delays from cross-correlation analysis. 
However we can establish quantitatively the evidence of contin- 
uum and, most important, of line variability. For this purpose we 
define the unbinned discrete structure function 



UDSF(Tu) 



\y(td-y(tj)\, 



(4) 



where v(f) represents any of the line or continuum light-curves 
considered, f, and tj are two observation epochs and r,y - f, - tj 
is the time delay. The (binned) structure function can be defined 
in bins of time delay, centered at t: 



Y^UDSFiTij) 



'.J 



(5) 



where the sum is extended to all the M values of UDSF belong- 
ing to a given bin of t. We adopt the average of UDSF, instead 
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days 

Fig. 4. Structure functions for PG 1634+706 in the observer's 
frame. Upper panels: emission lines, middle panel: shorter wave- 
length continuum Cshon, lower panel: longer wavelength contin- 
uum c,a„g, left: CIII1M1909 A, right: Mgll ^2798 A. 

of UDSF^, since it is more stable, i.e. less sensitive to deviant 
points; the factor in equation (4) makes S{t) equal to the 
standard deviation in the case of a Gaussian distribution (see|Di| 
IClemente et al.|[T996| l. Structure functions for lines and continua 




200 400 600 800 1000 200 400 600 800 1000 

days 

Fig. 5. Structure functions for PG 1247 +268. Upper panels: 
emission lines, middle panel: shorter wavelength continuum 
Cshort, lower panel: longer wavelength continuum ciong, left: 
CIV/il549 A, right: CIII] ^1909 A. 



of the two quasars PG 1634+706 and PG 1247+268 ai-e reported 
in Figures |4] and |5] respectively, as a function of the time delay 
in the observer's frame. The error bars reported in the figures 
represent simply the standard deviation of the UDSF in each bin 
and not the uncertainty on their average value. The latter would 
be 1/(M - 1) times the standard deviation, for M uncorrected 
UDSF values. 

In the case of PG 1634+706 the structure function analy- 
sis simply confirms what already appears from the light-curves, 
namely no significant variations of lines and continua are de- 
tected. The r.m.s. noise level is about 0.02-0.03 for continua and 
0.05 for lines. 

On the contrary, in the case of PG 1247+268, there is a clear 
increase of variability for long time lags. To establish the signifi- 
cance of this variability it is necessary to evaluate the probability 
of the null hypothesis that a value of S (t) is produced by pure 
noise. For this purpose we generated mock noise light-curves , 
n(ti), Nepo, where Nepo is the number of observing epochs, 
by extracting random numbers with a Gaussian distribution of 
standard deviation cr, which represents pure noise, under the as- 
sumption that the photometric noise is not correlated on the time 
scale of our minimum sampling interval (namely about 10 days). 
Then, the relevant structure function 5 mof/t(T) has been computed 
with the same binning adopted for real data. The simulation was 
iterated for 10^ times to generate, for each bin, the statistical dis- 
tribution of the S,„ock{T) values. 

The value adopted for cr has been estimated, conservatively, 
from the observed structure functions of figure |4] and |5] them- 
selves, as the value of S(t) in the first bin. This is, in fact, an 
overestimate of noise, since it includes short-time-scale varia- 
tions of the QSO, which however should be negligible on the 



basis of previous studies of (fainter) QSOs (Giveon et al. 1999 



Kaspi et al. 2000 1. The almost constant value of the structure 



functions in the case of PG 1634+706 confirms the adopted hy- 
pothesis that noise is not correlated on long time scale. 

For the CIII] line the rm.s. noise (~ 0.05) is larger than in 
the case of the CIV line (~ 0.02). This depends on the intrinsic 
faintness of the CIII] line. According to the above simulations 
the probability of the null hypothesis, for the second and fourth 
bins, is P(> S) ~ 4 x 10""*. The local maximum of S(t), for 
CIII] at T ~400 days, is due to an "oscillation" of the light-curve 
on a time scale comparable with the total time base of the ob- 
servations. This indicates that a longer time base is needed to 
properly average individual oscillations and to obtain a "station- 
ary" structure function. Further sampling, trying to avoid yearly 
periodicity due to periods of optimal observability, is needed to 
obtain the resolution required for time delay measurements. 

The structure function of CIV line shows a much more sig- 
nificant variability: the probability of the null hypothesis in the 
case of the last two bins, where S{t) is 0.13 and 0.21 respec- 
tively, is less than 10"^. The systematic decrease of CIV line 
luminosity and the relevant continua might seem to rely on the 
value of the light-curves at the first epoch. Thus we have recom- 
puted the S (t) after removing the first epoch. Indeed, in the case 
of continua the systematic decrease of S (t) is no longer evident. 
On the contrary, for the line it remains practically unchanged 
(except that the last bin is no longer present). Thus the fact that 
the probability of the null hypothesis is smaller than 10"^ does 
not rely on that particular point in the light-curve. 

The main goal of new echo-mapping campaigns is to obtain 
a direct measure of the central black hole mass and to establish 
whether the scaling of QSO broad line region size and masses re- 
main the same at higher quasar luminosity. The very detectabil- 
ity of line variation in the most luminous QSOs was, until re- 
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cently, an open question, since 3C273 with (/IL^(5100A) ^ 10' 
erg s~') was the brightest QSO with detected line variability 



(Kaspi et al. 2005 ). An evidence of line variability in brighter 
objects has been reported more recently by |Kaspi et al.| ( [2006| l. 
The present detection of line variability in one of the two quasars 
considered strengthens the notion that quasar emission lines re- 
spond to continuum changes in high luminosity QSOs like in 
lower luminosity ones and indicates that echo-mapping cam- 
paigns are feasible for ^L^(5100A) > 10"^^ 



erg s 



Our results are consistent with those of Kaspi et al. (2006i, 
who have monitored sources of comparable luminosities and 
redshifts. We have also detected variability, of comparable am- 
plitude, in the continuum and CIII] and CIV emission lines in 
one of our sources. Considering the shorter time baseline of 
our observations, the fact that variability has not been found in 
PG 1634-1-706 is consistent with the segments without variations 
which can be seen in the Kaspi et al. ( |2006| l light curves on this 
time scale. 

To estimate the further monitoring needed to measure line- 
continuum delays, we should take into account that, on aver- 
age, the amplitude of variability increases with the rest-frame 
time lag. However this dependency is weak: log S{Trest) ~ 
0.15 log Tres;, according to previous statistical studies on large 
quasar samples ( [Vanden Berk et al. 2004 ; De Vries et al. |2005| l. 
In the case of our quasars PG 1634+706 and PG 1247H-268, the 
observer's frame t„/,5 = Tresii^ +z) is dilated by a factor 2.3 and 3 
respectively. The size of the broad line region, based on Balmer 
line reverberation-mapping, is expected to scale with luminosity 
as Rblr with 7=0.5-0.7, depending on how L is measured 
and the Rblr vs. L fit is obtained (Kaspi et al. , 2000 2005"; Bentz' 
2006) l. A Balmer line BLR of about 3 



et al. 2006 Kaspi et al 



light years in the rest frame is expected for the relevant luminosi- 
ties. On the other the hand CIV line corresponds to a BLR size 
2 times smaller than Balmer lines (Peterson & Wandel 2000). 
Considering that in general a time baseline of about twice the 
light-crossing time is required in order to obtain a reliable lag 
from the light curves, we estimate that reverberation mapping 
should be feasible after 5-6 more years of monitoring. 
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4. Summary 

- We have initiated a campaign for the monitoring of 4 high 
luminosity QSOs and we present the result for two of them 
on the basis of 3.3 years of observation with the 1.82 m tele- 
scope of the Asiago Observatory. 

- We discuss the data reduction, the procedures adopted and 
the level of accuracy attained in relative spectrophotometric 
variability measurements. 

- We perform a structure function analysis of the light-curve 
and numerical simulations to establish the confidence level 
of line variability detection. 

- We detect line and continuum variability in one of the two 
QSOs, PG 1247+268 of .1L^(5100A) = 10^^ erg s"', with a 
probability of less than 10 of the null hypothesis that the 
observed structure function is produced by pure noise. 

- This detection supports the notion that emission lines re- 
spond to continuum variations as in substantially less lumi- 
nous QSOs. 

The results encourage the prosecution of the campaign which 
should provide time delay, and black hole mass, estimates in 5-6 
years. 
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